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1 Introduction 



In a homogeneous, isotropic, expanding universe, the apparent luminosity of a standard 
candle is a monotonically decreasing function of the time of flight of emitted photons 
if the universe is open. This is also true in spherical universes if the time of flight is 
small enough with respect to the radius divided by the speed of light. In principle the 
(apparent) luminosity £{t) as a function of time can be used to measure the scale factor 
a{t). In reality, arriving photons do not tell us their time of flight, but only their spectral 
deformation, z := z^tabie/t'observed — 1- In an expanding universe z is positive, 'red shift'. If 
we pretend to know the scale factor a{t) we can compute the luminosity i{z) and confront 
it to the Hubble diagram. The transform a{t) i{z) reminds us of the Fourier transform 
and of course we are interested in the inverse transform i{z) — > a{t). Therefore we must 
ask three questions: What is the domain of deflnition of the initial transform, what is its 
image and is the transform injective? Should the measured luminosity i{z) be 'far away' 
from the image our working hypotheses are put to test. 



We assume the kinematics of general relativity: (i) The gravitational field is coded in 

a time-space metric of signature H , the configuration space is the set of all such 

metrics, (ii) Massive and massless, pointlike test particles, subject only to gravity, follow 
timelike and lightlike geodesies. (Hi) Pointhke clocks, e.g. atomic clocks, are necessarily 
massive. They move on timelike curves (not necessarily geodesies) and indicate proper 
time r. 

This kinematics is covariant under general coordinate transformations. In 1983 the 
meter was officially abandoned as fundamental unit in favor of an absolute speed of light 
and we would like to stress that at least since then any non-covariant kinematics is void 
of any physical meaning. 

We also assume the hypotheses of spatially maximally symmetric cosmology: (iv) The 
metric is Robertson- Walker, 



The scale factor a{t) is a strictly positive function of time, k — —1, 0, 1 for the pseudo- 
sphere, the Euchdean space and the sphere. We take the coordinates x dimensionless and 
call them 'co-moving position', while the scale factor is measured in meters. For closed 
universes, k = 1, we must restrict the spatial coordinates to the unit ball, < 1, they 
describe the northern hemisphere. For open universes, x varies in M^. (v) The test parti- 
cles are (super clusters of) galaxies and photons. The former follow the timehke geodesies 
t — T, X — constant. Note that due to the high degree of symmetry the proper time is 
universal for all these timelike geodesies and is taken as time coordinate. The photons 
are emitted from the a galaxy at time t and arrive at our position at time to, today. 



2 The hypotheses 




(1) 
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The symmetry hypothesis and the choice of the coordinates {t, x) have reduced the 
configuration space to the set of positive functions a{t) on the real fine and to the para- 
meter k = 0, ±1. 



3 The Hubble diagram 

^From these hypotheses we can compute [T] the (apparent) luminosity i{t) of a standard 
candle in Watt/m^ ClS db function of emission time t 

m = -± (2) 

where L is the absolute luminosity of the standard candle in Watt, := a(to) is the scale 
factor today, 

cdt 

Xit) := / -J- (3) 
Jt a{t) 

is the dimensionless co-moving geodesic distance covered by the photon and 

f sinx, k = 1 
six) ■.= { X, k = . (4) 




The luminosity has a singularity at t = Iq, ' short distance divergency', for any value of 
k and any scale factor. For closed universes, = 1 we may have additional singularities 
for X = vr, 'antipode divergencies', for x = 27r we are back at a second short distance 
divergency and so forth if x continues to grow. Of course there might be a horizon, i.e. 
an upper bound for x, masking all or some of the singularities. 

/^From the above hypotheses we can also compute the spectral deformation as function 
of emission time, 

m = - 1- (5) 

a[t) 

The theoretical Hubble diagram is the parametric plot in the z — £ plane as t varies. 

If we suppose that the scale factor is strictly increasing with a := da/dt > then z is 
positive, 'red shift', and we can invert the function z{t). By abuse of notations we write 
t{z) for its inverse. Then the Hubble diagram is a function i{t{z)) that still by abuse is 
written £{z), 



where 



X{z) := X{t{z)) = - f 4k 
ao Jq H[z) 
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and H{z) := a{t{z)) / a{t{z)) is the Hubble rate. 

The short distance divergency now is at ^ = and easy to get rid of: let us define the 
regularized luminosity 

f{z) := z^£{z) = , . (8) 

Indeed, 

where the prime is differentiating with respect to z, q{z) is the deceleration parameter, 

g(z) = -^(t(z)), (10) 
aa 

and Ho := i/(0), go = g(0). 



4 The transform a{t) — > £{z 

Let us first try to describe the image, that is all luminosity functions i{z) which can be 
obtained from strictly increasing scale factors a{t) with d > and with k = 0, ±1. 

We already know that i{z) comes with the short distance divergency at 2; = which 
is such that z'^i{z) =: f{z) is regular there. For open universes there are no other sin- 
gularities. Indeed, {z + l)^i{z) =: g{z) is a decreasing function. For closed universes on 
the other hand, g{z) goes through a minimum as the photons pass the equator, x = 7r/2. 
From there on the luminosity increases again and goes to the antipode divergency. It 
might of course happen that the equator is masked by the horizon in which case g{z) 
remains decreasing for ever, even though the universe is closed. 

Let us now ask whether the transform is injective in the domain of increasing scale 
factors. 

If we pretend to know k the answer is affirmative for open universes. Indeed, solving 
equation (jS)) for s{x) and differentiating with respect to z yields 

«o . , iz + l),/A7TalL-^iz + iyf{z)-kz'^ _ z{z + l)^AT:alL-^g{z) - k 

Therefore we can reconstruct the Hubble rate from the luminosity. 

Integrating the Hubble rate with respect to t gives us the scale factor with the ambi- 
guity of the initial condition oq. But for flat universes this initial condition is unphysical, 
by a coordinate transformation of x, we can set oq = 1 m. This is different for curved 
universes where is related to a local observable, curvature. In the closed case, ao is 
also related to a global observable, the radius of the universe today. However, unless g{z) 
already exhibits an increase, only a lower bound on the radius today can be reconstructed 
from the luminosity. 



ao>\ -A ' f7min:=min{(z + l)2£(2;)}. (12) 



4 



Note that this lower bound does not depend on the absolute luminosity L. 

If we admit that we do not know k and if the luminosity function i{z) satisfies: (i) 
z'^i{z) is regular at z = 0, {ii) {z + lYi{z) is decreasing, then there are three positive 
functions a±(t) and a(t), such that the universes with scale factor a_(t), k = —1, with 
scale factor a{t), k = and with scale factors a+{t), k = 1 have the same luminosity 
function i{z). These three scale factors satisfy 

a_(to) sinhx- = a(to) X = «+(io) sinx+. (13) 

Note that in the flat case the 'initial condition' a(to) carries no information while in the 
closed case a+(to) must satisfy the inequality (fT^ . 

Example (constant deceleration parameter): 
Take the scale factor, 

a{t) := aoipHoty/P, p > 0, < t < to = 4r • (^4) 

Then the Hubble rate is H = Hq{z + 1)^ and the deceleration parameter is constant, 
q = p ~ 1. The co-moving distance is 

X= ^ln{z + l) = -^ln{Hot), p = 1, (15) 



X 



{p - l)aoi/c 



{l-{z + lf-^) 



{p - l)aoifo 



:i-{pH,t) 



i-i/p\ 



p^l. (16) 



For p > I, there is a horizon at x = c/ {{p — l)ao-ffo)- In particular for p = 2, k = the 
regularized luminosity f{z) is constant. 

Example (constant regularized luminosity): 

Suppose we have measured a constant regularized luminosity f{z) = LHq/{Atic^). Then 
the Hubble rate is 



H={z + l)\ H^{z + 1)2 - kc'^zyal 



(17) 



The three solutions of these three differential equations, k = ±1, 0, in terms of the scale 
factors are obtained from equations (fT^ and ()14|) : 



a{t) = ao{2HQty/^\ll- k 



1 - {2Hoty/^ 
aoHo/c 



0<t<t 



1 







2Ha 



To alleviate notations we have suppressed the subscripts from a{t) and oq. In the closed 
case, oo must satisfy the inequality (fT^: 



ao > 



c z„ 



Hq Z^ax + 1 



(19) 
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and the inflection point of g{z) is hidden behind the horizon at x = arcsin[c/(ao-ffo)]- 



Example (constant Hubble rate): 

To have an example without horizon consider the scale factor 

a(t) = aoexp[Ho(t - to)]. 



(20) 



This is a limiting case of the first example with p 1. It has constant Hubble rate 
and constant deceleration parameter, H{z) = Hq, q{z) = —1, and x = cz/{aoHo). In the 
closed case, g{z) has an infinite number of inflection points alternating with short distance 
divergencies. The inflection points are all minima, the first is located at z^fi = vrao-f^o/ (2c). 



Example (closed universe): 

Suppose we have measured the luminosity: 



LHl z + 2 
Attc^ z'^ 



(21) 



It is strictly decreasing and has the correct short distance divergency. Its g{z) has one 
and only one inflection point at Zi^a = (1 + vT7)/2 ~ 2.56 and suggests a closed universe 
with ao ~ 0.337c/ifo- 

Counter-example (wiggling g{z)): 

Suppose we have measured the luminosity, figure 1: 



i{z) 



LHl (sinz/;z)^ + 0.1 
47rc^ z"^ 



(22) 



Again it is strictly decreasing and has the correct short distance divergency. Now g{z) has 
a maximum, figure 2. Therefore no Robertson- Walker universe, neither open nor closed, 
exists with this luminosity. 



Figure 1: The monotonic luminosity 



The last two examples illustrate that our constraint of monotonic g{z) is stronger than 
the constraint of monotonic luminosity (!-{z). 
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Figure 2: Its wiggling g{z) 



5 Non-monotonic scale factors 

If the scale factor is strictly decreasing we get similar results with a negative spectral 
deformation: — 1 < 2; < 0, 'blue shift'. 

One might think that one can produce non-monotonic functions g{z) by starting from 
non-monotonic scale factors a{t). This is not true. In fact any non-monotonic scale 
factor produces multivalued luminosities in terms of the spectral deformation z. The first 
example is of course the constant scale factor with no spectral deformation, z = 0, but 
varying luminosity. A more generic example is, see figure 3: 




Figure 3: The non-monotonic scale factor (j^Hj) 



a{t) = lgr + a, g,a>0. (23) 

For positive t^, we have a maximal redshift of z^i,^ = ao/« — 1. We have z = for t = ±to 
and a blue shift for t < —tQ-. 

- 1 < z^,^ = 1< 0. (24) 

^gt'' + a 
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The dimensionless co-moving distance is 



X{t) = c\ — 
ga 



arctan 



arctan 



At 
2a 



(25) 



and the relation between emission time t and spectral deformation z is 



t(z) 



± 



ao 



z+1 



a. 



(26) 



Note that the Hubble rate vanishes at the inflection time t = 0, z = 2;max while the 
deceleration parameter diverges there. Note also that the regularized luminosity vanishes 
at t = —to, z = 0. As t tends to minus infinity z tends to —1 and the luminosity, 
regularized or not, diverges, the 'ultra-violet divergency'. In the closed case, we have 
in addition antipode and short distance divergencies. Figure 4 shows only one antipode 
singularity. 



/ 



^ ~ ^ LHQ/(4nc) 




t - 



Figure 4: The regularized luminosity of the non-monotonic scale factor 



6 Data and conclusion 

We use the 'Gold' sample data compiled by Riess et al. (2004) [3 , with 157 SN's including 
a few at 2; > 1.3 from the Hubble Space Telescope (HST GOODS ACS Treasury survey). 
For convenience we normalize the luminosity to the maximum absolute SN luminosity 
estimated by Jha et al. (1999) gj, Saha et al.(2001) and Gibson & Stetson (2001) 
[HI, L = L^ax = (1 ± 0.1) ■ 10^^ W. The Hubble rate today Hq is taken as (70 ± 5) km 
s"^ Mpc-i [ZIE!. 
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Figure 5: The regularized luminosity as measured today [S], with a binning of 0.02 in red 
shift. The full line at low red shift corresponds to a second order polynomial extrapolation 
fit. 



The regularized luminosity allows to extract LHq from the Hubble diagram with small 
red shift, figure 5. The value /(O) is extracted by a second order polynomial extrapolation 
fit on the SN data up to a red shift of 0.1. By construction, the fitted value /(O) is equal 
to LHq/{4:7ic^) = (4.3 ± 0.4) ■ 10^^ W Mpc~^ where the error is only coming from the fit 
itself. 



6.1 Lower bound on the radius of the universe today 

The first of equations ^ and equation (fT^ give the lower bound on oq as 




ao>\^—. (27) 



The minimal value of g{z) is obtained from the SN recorded at z = 1.3 and is equal to 
9min = (3.83 ± 0.72) • 10^6 W Mpc^^. At 95% confidence level (CL) we have the upper 
bound on the minimal value of g{z): Qnim < 5.6 ■ 10^^ W Mpc""^. Combining the errors 
from /(O) and g^in gives the lower bound on gq at 95% CL: 

ao > 0.88 c/Ho ~ 3.8 Gpc ~ 1.2 ■ 10^^ m ~ 1.3 ■ 10^° light years. (28) 
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This lower bound on ao translates into a lower bound on the curvature density, 

Qj^---k(-^\ > -1.29 at 95% CL. (29) 
V (^oHq J 

Note that this limit is independent of the values of the absolute luminosity L and of the 
Hubble rate today. It is also independent of any dynamical hypothesis. 

Let us compare this bound with the one obtained from the SN data but now adding 
the dynamics of the ACDM cosmology fitting the matter density Qm, the cosmological 
constant density Q\ and the nuisance normalisation parameter and without any other 
input constraint: > at 95 % CL. 

6.2 Wiggles and non-monotonicity 

We must now ask the question whether the data is compatible with a monotonic luminosity 
i{z). We will also ask the finer question whether g{z) is monotonic. 

To detect non-monotonicity in the SN data set we assume that the luminosity i{z) 
and its g{z) can be described by monotonic functions to which we add a simple Gaussian: 

L{z) = i{z,LHlp) +a^exp-{z - z^)y{26zl) (30) 

and 

G{z) = giz,LHlp) + a^exp-iz - z^Y/i25zl). (31) 

The monotonic functions i{z, LHq,p) and g{z, LHq,p) are derived from a power law 
parameterization of the scale factor a{t) := ao{pHoty^P (constant deceleration parameter), 
l/5<p<l, fc = 0. It yields the monotonic luminosity 

and with p = 0.69 it describes well the ACDM cosmology up to a redshift of 1.8 [Qj. Its 
acceleration is positive with q = —0.31. 

Our wiggle detection procedure consists of scanning the plane in wiggle position z^ 
and wiggle width in steps of 0.01 in both directions. In each point (z^,^^;^) of this 
plane we fit the normalization LH^, the power p and the wiggle amplitude a^- 

Warning: if the wiggle amplitude is smaller than a critical amplitude the mod- 
ified functions (|30p and (|31|) will still be monotonic. 

We would claim that i{z) and a fortiori the luminosity g{z) is not monotonic if the ratio 
between the fitted wiggle amplitude and the associated error is greater than 5 {5a level 
detection) and if the wiggle amplitude is greater than the critical one. The sensitivity of 
the method is computed by Monte Carlo simulation. The same SN sample than the Riess 
data set with the same statistical power is simulated assuming the ACDM cosmology and 
a wiggle of positive or negative amplitude is added to i{z) and g{z) for each point in the 
{zw,Szw) plane. We apply the wiggle detection procedure on each simulation, restricted 



10 



to a small grid of points around the simulated one to speed up the processing. The 
significance on the wiggle amplitude ( fitted) I /o"a„) is computed in each point and the 
smaller value from positive or negative wiggle amplitude is retained. The sensitivity is 
computed at a 1o level corresponding to a 95% confidence level exclusion limit on the 
wiggle magnitude defined by 

a^{z^) = ±(l0-"»/2-5 - iy(zj. (33) 

Figure 6 shows the significance of the wiggle fit performed on the Riess data sample 
(color contours) for the luminosity i{z) and g{z) with Zw varying from 0.01 to 1.8 and 
Szm from 0.01 to 2 in steps of 0.01 in both directions. The maximum significance for both 
i{z) and g{z) is 2.4 for a wiggle at the position z^ = 0.45, with width of 0.07. The wiggle 
magnitude is = 0.16 for the luminosity i{z) and = 0.25 for g{z). The dashed lines 
indicate the location of the critical wiggle magnitude of a positive wiggle that breaks 
the monotonicity. 

No wiggle greater than this value is observed and we conclude that no wiggle is detected 
at a 5a level using the actual SN data set. On the same figures, the sensitivity for different 
values of the wiggle magnitude is shown (plain line) . Wiggles of magnitude greater than 
2 are excluded at 95% CL up to a redshift of 1.6. Up to a redshift of 1, the 95% CL 
exclusion limit on the wiggle magnitude is 0.6. These two magnitudes are below the 
critical magnitudes and therefore these wiggles do not upset the monotonicity. 
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Figure 6: Color contours: Significance of wiggle detection (vertical colour scale) as func- 
tion of the wiggle position Zy^ and width 5^;^, logarithmic scale. Full lines: Expected 2a 
sensitivity of the wiggle detection as a function of wiggle magnitude m^. The dashed lines 
indicate the location of the critical wiggle magnitude of a positive wiggle that breaks 
the monotonicity. The upper panel is the exclusion plot for a non-monotonic luminosity 
^{z), the lower panel for non-monotonic g{z). 



